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Abstract 

Using a recently developed formalism of quantization of radiation in the presence of 
absorbing dielectric bodies, the problem of photon tunneling through absorbing bar- 
riers is studied. The multilayer barriers are described in terms of multistep complex 
permittivities in the frequency domain which satisfy the Kramers-Kronig relations. 
From the resulting input-output relations it is shown that losses in the layers may 
considerably change the photon tunneling times observed in two-photon interference 
experiments. It is further shown that for sufficiently large numbers of layers interfer- 
ence fringes are observed that cannot be related to a single traversal time. 



PACS number(s): 42.50.Ct, 73.40. Gk, 42.79.-e 



1 Introduction 



Stimulated by recent experiments H ||, the problem of photon tunneling through 
multilayer dielectric barriers has been of increasing interest. In order to answer the 
question of what is the time that is spent by a photon inside such a barrier, the effects 
of dispersion and absorption should be considered very carefully. The calculations 
that have been performed so far are based on real refractive indices of the layers 
El 0, ||, so that a number of questions, such as the influence of absorption on 
the measured traversal times 0, have been open. It is well known that in frequency 
intervals where the dielectric layers are nearly transparent the action of each layer and 
hence that of a multilayer barrier can be described in terms of unitary transformations 
that relate the operators of the outgoing fields to those of the incoming fields (see, 
e.g., Ref. [§]). These input-output relations and the underlying quantization scheme 
(see, e.g., Refs. [[7], ||, |J) of course fail, when the effects of absorption cannot be 
disregarded. 

Various approaches to the problem of quantization of radiation in the presence of 
absorbing dielectric bodies have been developed [jlUJ [TJ, 0, [T3|, [14], [15|, [19, O, [IB], |1§ 



20| , p| . In the present paper we use the quantization scheme given in Refs. [^0] |2l|| . 
It is based on a Green function expansion of the operator of the (transverse) vector 
potential and applies to radiation in both homogeneous and inhomogeneous dielectric 
matter. In this approach, the matter is described in terms of a complex permittivity (in 
the frequency domain), without using a particular microscopic model for the matter. 
The only condition is that the permittivity satisfies the Kramers-Kronig relations, 
because of causality. Applying the method to the calculation of input-output relations 
for radiation at absorbing multilayer dielectric barriers, we can systematically study 
the effects of dispersion and absorption on the propagation of single-photon wave 
packets through such barriers. 

The paper is organized as follows. In Sec. ^]the quantization scheme is outlined. In 
Sec. |3| the scheme is applied to radiation falling on multilayer dielectric barriers. The 
input-output relations derived are used in Sec. f| in order to calculate barrier traversal 
times measurable in two-photon interference experiments. Finally, a summary is given 
in Sec. ||. 



2 Field quantization 

Let us consider linearly polarized light propagating in a; direction in an inhomogeneous 
linear dielectric medium. Introducing the (transverse) vector potential 

A(x,t)= du e~ iuJt A(x, u) + c.c. (1) 



1 



(A = A y ), the classical (phenomenological) Maxwell equations yield 



d 2 uj 2 
— + — e(x,uj) 
ox c 



A{x,uj) = 0, 



where 



e(x, uj) = e r (x, uj) + i e ; (x, u) 



(2) 



(3) 



is the complex permittivity which for inhomogeneous media, such as multilayer di- 
electric barriers, varies with x. Clearly, when the permittivity is complex, then Eq. (|2|) 
cannot be valid as an operator equation in quantum theory. On the other hand, it is 
well known that propagation of light in absorbing matter is unavoidably accompanied 
by noise. In a quantized version of Eq. (|2|) this noise source can be taken into account 
by introducing an operator noise current j n , so that 



d 2 



UJ 



+ — e{x,uj) 



dx 2 c 2 



A(x,uj) =j n (x,u), 



where 



UJ 



h 



£i(x,Uj) f(x,u). 
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Here, f(x,uj) is a bosonic basic field, 

f(x,uj)J\x',uj')\ = 5(x - x')5(uj - uj'), 

f(x,u)J(x',uj') = P(x,uj), }\x\uj') 



0. 



(4) 



(5) 



(6) 
(7) 



and A is a normalization area in the yz plane. Equation (|]) is now an equation for 
the the (Schrodinger) operator A(x, uj). The solution can be represented as 



A(x,uj) = / dx' G(x, x' , uj)j n (x', uj), 



(8) 



where G(x,x',uj) is the classical Green function that satisfies the equation 



d 2 uj 2 
— - + —e{x,uj) 
ox c z 



G(x, x', uj) = 5(x — x') 



(9) 



and vanishes in the limit when x— >±oo. The quantization scheme ensures that the 
well-known equal-time commutation relation 



A(x,t),E(x',t) 



ifi 



At { 



5(x — x) 



(10) 



is preserved PT 
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3 Input— output relations 

Let us consider a dielectric barrier consisting of (N — 2) layers (N>3), 



N 



e(x,uo) = J2 X j( x ) € j 



to 



(11) 



3=1 



where ej(u) is the permittivity of the jth layer, and 



Xj(x) 



1 if Xj-i < X < r 
otherwise 



(12) 



(x — > — oo, xat^oo). From Eqs. (P, (||), and (||), with e{x,u) from Eq. (|Tl[), the 
operator A(x,u>) can be rewritten as |2l|] 



K J fr{ Jo \ ±ircue ej(uj)A 



+ H.c. 



(13) 



where the quasi-mode operators a J + (x,u;) and aj_(x,c<j), which are associated with 
the (damped) waves propagating to the right and left, respectively, are related to the 
bosonic basic field as 



a j± {x,to) = cij± [x ,to) e 
1 / ~~Zo 



=F7j (lo)uj(x— x')/c 



I 



±-J2ij(u)- / dye 



C J x 



f(y, w) 



(14) 



(xj-i<x,x' <Xj). Here the notation J ej{to) = rij{to) = (3j(to)+i / yj(to) is introduced. 

Using the commutation relation (^|), from Eqs. (|TJ) we find that the quasi-mode 
operators a 1+ (x,to), x<x lt and a N _(x,to), x>x N , of the incoming radiation from 
the left and the right, respectively, satisfy the commutation relations 



a 1+ {x, to), a\ + {x', to')] = S (to - to') e -'nM«l*-* , IA ) 
a N _{x,to),a ] N _{x',to')] = 5 {to - to') e^^MI/^ 



a 1+ (x,w),alf_(x',u') 



(15) 
(16) 

0. (17) 

Note that Eqs. ( |i~5l) and ( jlBf) agree with the commutation relations valid for the 
corresponding bulk materials. For vanishing absorption (7i(at)(^) — ► 0) the operators 
a 1+ and a N _ are independent of a; and ordinary free-field bosonic operators | 2l|| . 

The output operators ai_(x,uo), x < xx, and d]y + (x,to), x>x^, can then be 
calculated step by step starting from a single-slab plate (iV = 3). Using Eqs. ((13|) and 
( P^| ) and taking into consideration that the vector potential must be continuously 



differentiable at the interfaces, after some lengthy calculation we find that the output 
operators can be expressed in terms of the input operators and bosonic operator noise 
sources g±(u) associated with the losses in the barrier 0, pBfl , 



ai_(a; 1 ,a;) 
a>N+( x N-i, 



TM 



a 1+ {xi,u) 




g±{w),g±{u') =5(w-iJ) 



(18) 
(19) 



(the input and noise operators are commuting quantities). Here the characteristic 
transformation matrix T(a>) describes the effects of transmission and reflection of 
the input fields [^2|, whereas the losses inside the barrier give rise to an absorption 
matrix A(cu). Explicit expressions for the matrices T(lu) and A(cu) and the noise 
operators g±{uS) [as linear functional of the field f(x,u) inside the barrier] are given 
in Ref. IE3I. 



The input-output relations (0) [together with Eq. (|i~4|)1 lead to commutation re- 
lations for the output operators ai-(x, to), x<X\, and a N+ (x, u>), x>xm, that differ, 
in general, from those in Eqs. ([15]) - ([17]) for the input operators. The difference de- 
creases with increasing distances from the barrier. In particular, it can be disregarded 
when the sourrounding medium can be regarded as beeing lossless and the input and 
output operators are ordinary bosonic operators. In this case the relations 



|TnH| 2 + |T 12 H| 2 + \A xl (u})\ 2 + \A 12 (u)\ 
= \T 21 (lu)\ 2 + \T 22 (u)\ 2 + \A 21 (cu 
T,J,n\T* -i- T^(i,i\T*d,i\ -i- A, J,, A/ 



L21 HI 2 + |a 22 hi 2 = i, 

T x1 {u)T; x {u) + T 12 ( W )r* 2 ( W ) + A u (uj)A* 21 (uj) + A 21 (lu)A* 22 (lu) 



(20) 
(21) 



can be shown to be valid, which ensure preservation of the bosonic commutation 
relations. 



4 Photon tunneling 

To study the influence of dispersion and absorption on photon tunneling through 
multilayer dielectric barriers, let us consider a two-photon experiment of the type 
described in Ref. |l]] (Fig. 1). Pairs of down-conversion photons are directed by mirrors 
to impinge on the surface of a 50%:50% beam splitter and the output coincidences 
are measured. One photon (I) of each pair travels through air, while the conjugate 
photon (II) passes a barrier. The coincidences attain a minimum when the two 
photons' wavepackets overlap perfectly at the beam splitter. This can be achieved by 
translating an appropriately chosen prism in one arm of the interferometer in order 
to compensate for the delay owing to the barrier. 

Let us assume that the barrier is in the ground state and the two correlated 
photons are prepared in the state 

1^)=/ dQa{Q) dufi^fiQ-^ali^aliQ-u)^), (22) 
Jo Jo 
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where a(Q) and f{uj) are the bandwidth functions of the laser and down-conversion 
photons, respecteively, f(tu) being centered at Vt/2. From photodetection theory it 
is well known (see, e.g., |25[]) that the overall coincidences R can be obtained from 
the time-integrated normally ordered intensity correlation function, 

R = eJdt 1 J dt 2 (^-)(t 1 )^ H fe)^ (+) (ii)^ (+) fe)), (23) 

where E^(ti) and E^{t 2 ) are the fields at the detectors in the two output channels 
of the beam splitter (£, detection efficiency). Applying the input-output relations (|18|) 
and using Eq. (^), after some lengthy but straightforward calculation we find that 

/•CO 

R = 2tt 2 AT 4 / dtta 2 (n) F(n), (24) 
Jo 

F(n) = [ n duj f 2 (uj)f 2 (n-uj) uj(^Uu)T* (iUu)\T 12 (n-uj)-e- 2ins e 4iuJS T 12 (uj)}, 
Jo 1 J 

(25) 

where s is the translation length of the prism (cf. Fig. 1), and the abbreviation M — 



£h/(4:7rce A) has been introduced. 
The translation length s= so that corresponds to the minumum of R(s) is usually 
used in order to distinguish between superluminal (positive values of so) and sublu- 
minal behaviour (negative values of s ) of the photon passing through the barrier. In 
the numerical calculations we have considered H(LH) fe structured plates (H, titanium 
dioxide; L, fused silica) of A/4-layers, which are of the type described in Ref. H|. We 
have calculated the function T(uj) applying the algorithm given in Ref. []23[]. Due to 
the lack of reliable data we have assumed that in the (relevant) frequency interval 
the complex refractive indices are approximately independent of frequency, so that 
all the dependences on frequency effectively result from the geometry of the barrier. 
For the sake of simplicity we have assumed that the line shape function of the ex- 
citing laser, a(Q), is sufficiently small, so that F(VL) f»F(u;o) in Eq. (p4|), where lu 
is the centre frequency (uj = 5.37 x 10 15 s _1 ). Introducing the single-photon pulse 
shape function f(t) = (2ir)~ 1 ^ 2 J dcjexpf— iut}f(u), we have performed calculations 
for both Gaussian pulses f{t) ocexp[iu; t/2 — (t/t ) 2 } and time-limited non-Gaussian 
pulses f(t) oc exp{icu t/2 - [1 - [t/(2t )] 2 ]- 1 } if |t| < 2t and f(t) = elsewhere, 
where t = 20fs in either case. 

The values of Ar = 2so/c that are shown in Fig. 2 are valid for both Gaussian and 
time-limited non-Gaussian pulse shape functions. The values are positive and indicate 
superluminal behaviour of the photon at the barrier, the characteristic tunneling time 
being given by r t = l/c — At (I, thickness of the barrier). From the figure we see that 
the "lead" of the photon, Ar, increases with the number of layers of the barrier, N 
= 2k + 1, and tends to a linear function of N, the slope of which sensitively depends 
on the losses in the barrier. Disregarding the losses, the slope in the linear regime 
is simply given by the inverse velocity of light in vacuum, which indicates that r t is 
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independent of N. The effect of losses is seen to decrease the slope which implies 
that r t is increased. 

The interval of N in which Ar (linearly) increases with N must of course be 
limited by an upper boundary, which may substantially change with the pulse shape 
function of the photon at the entrance plane. For the system under consideration the 
increase of Ar with N ends when N^35 (lossless case) or iV^41 (lossy case) for the 
time-limited non-Gaussian pulse, whereas for the Gaussian pulse the boundary value 
of N is substantially increased. It should be pointed out that the observed increase 
of Ar with N is not in contradiction to causality. The effect can simply be explained 
by a shift of the pulse maximum towards earlier times owing to pulse reshaping in the 
barrier, where the earliest time is given by the time at which the pulse starts from 
zero. Since in the case of a Gaussian pulse the pulse maximum can be shifted to 
earlier times than in the case of a time-limited pulse, in the former case the upper 
boundary of N is higher than in the latter case. 

For sufficiently large iV the introduction of the time r (and the traversal time r t ) 
makes little sense. In Fig. 3 the transmittance of the barrier, T 12 (u), is plotted for 
relatively low (11) and high (41) numbers of layers. Since the spectral line shape 
function of the outgoing photon, f(uS) oc f(u;)T 12 (u}), sensitively depends on the 
two competing quantities f(uo) and T 12 (uj) (Figs. 4, 5), it can essentially be different 
from that of the incoming photon when the number of layers is large enough (compare 
Fig. 4 with Fig. 5). The behaviour in the time domain is illustrated in Figs. 6 and 
7 in which the intensity of the outgoing photon, I(t) = (E^E^), is plotted. In 
particular, for sufficiently large iV the incoming and outgoing photons' wavepackets 
lose all resemblance to each other (cf. Figs. 5 and 7). In this case the measured 
coincidences are expected to be more or less complicated functions of the translation 
length, the structure of which does not allow one to define uniquely a traversal time. 

Figures 4-7 refer to the case when the pulse of the incoming photon is time- 
limited. Compared to such a pulse, the spectral line shape function f(u) of a Gaussian 
pulse is more smoothed and its wings decrease substantially faster. Hence, the trans- 
formed line shape function f(u)T 12 (uj) of a Gaussian pulse reflects the frequency 
response of the transmittance of the barrier, T 12 (uj), less sensitively than that of a 
time-limited non-Gaussian pulse. This different behaviour explains the above men- 
tioned difference in the boundary values of N. 

In Figs. 8 and 9 the coincidences are shown as a function of the translation length 
for the chosen time-limited pulse shape and various numbers of layers. We clearly see 
that when the value of iV exceeds an upper boundary value, then the function R(s) 
loses the simple structure that can typically described by a well-defined minimum. 
It should be noted that, compared to lossless barries, frequency-selective absorption 
shifts the boundary towards higher values. With increasing N interference fringes 
are observed, which correspond to the various possibilities of (partial) overlapping of 
the undisturbed and the multi-peaked outgoing photons' wavepackets at the beam 
splitter. Clearly, each minimum introduces its own characteristic time, and a unique 
tunneling time can be hardly derived in this way. 
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5 Summary and Conclusions 



On the basis of a Green function approach to the problem of quantization of ra- 
diation in inhomogeneous, dispersive and absorptive linear dielectrics we have de- 
rived quantum optical input-output relations for optical fields at multilayer dielectric 
plates, which can be regarded as generalizations of the well-known concepts of unitary 
transformations that apply to non-absorbing matter. Applying the theory to photon 
tunneling through absorbing barriers, we have shown that relatively small imaginary 
parts of the refractive indices of the layers can already give rise to observable effects 
in two-photon interference experiments as performed recently. 

The results reveal that only up to an upper boundary for the number of layers 
the measured coincidences can be used for extracting from them a characteristic time 
that may be regarded as traversal time. The boundary value sensitively depends on 
the spectral line shape function of the photon at the barrier and the dependence on 
frequency of the transmittance of the barrier, which can be substantially different for 
absorbing and non-absorbing barriers. It is worth noting that for sufficiently large 
numbers of layers the photon's wavepacket can be distorted in the barrier to such an 
extend that the observed coincidences show a number of interference fringes which 
correspond to different time constants. 
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Figure 1: Scheme of the two-photon interference experiment QX|, |] for determining 
photon traversal times through multilayer dielectric barriers (L, laser; P, prism; DB, 
dielectric barrier; BS, beam splitter; PDi, PD 2 , photodetectors). 
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Figure 2: The temporal "lead" At = 2s /c that corresponds to the position s of 
the minimum of R(s) is shown as a function of the number of layers, N = 2k + 1, 
for a H(LH) fc structured plate of A/4-layers of the type described in Ref. [|l]; curve 
(1): lossless barrier (n T io 2 = 2-22, n Si02 = 1.41), curve (2): absorbing barrier (n Ti02 
= 2.22, n sl0a = 1.41 + 0.0372 1 
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Figure 3: The square of the absolute value of the transmittance of a multilayer 
absorbing barrier, |T 12 (co , )| 2 , is shown for AT = 11 layers [curve (1)] and iV = 41 layers 
[curve (2)]. The data are the same as in Fig. |[ 
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Figure 4: The (normalized) spectral line shape function f(uj) (full line) of a photon 
after having passed through an absorbing barrier consisting of iV = ll layers (data as 
in Fig. 0) is shown. For comparison, the line shape function of the incoming pulse 
that is assumed to be limited in time (2t = 40fs) is also shown (dashed line). 
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Figure 5: The (normalized) spectral line shape function f(u>) (full line) of a photon 
after having passed through an absorbing barrier consisting of iV = 41 layers (data as 
in Fig. 0) is shown. For comparison, the line shape function of the incoming pulse 
that is assumed to be limited in time (2t = 40fs) is also shown (dashed line). 
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Figure 6: The (normalized) intensity I(t) (full line) of a photon after having passed 
through an absorbing barrier consisting of N = 11 layers (data as in Fig. ^]) is shown. 
For comparison, the intensity of the incoming pulse that is assumed to be limited in 
time (2t = 40fs) is also shown (dashed line). 
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|I(t)| (normalized) 
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Figure 7: The (normalized) intensity I(t) (full line) of a photon after having passed 
through an absorbing barrier consisting of iV = 41 layers (data as in Fig. ^]) is shown. 
For comparison, the intensity of the incoming pulse that is assumed to be limited in 
time (2t = 40fs) is also shown (dashed line). 
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R (normalized) 




Figure 8: The (normalized) coincidences R(s) are shown in dependence on the 
translation length s for a time-limited pulse of the incoming photon (2t = 40fs) and 
various numbers of the layers of a lossless barrier: AT = 11 (dotted-dashed line), iV = 
35 (full line), N = A1 (dashed line). The data of the lossless barrier are the same as 
in Fig. |. 
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R (normalized) 




Figure 9: The (normalized) coincidences R(s) are shown in dependence on the 
translation length s for a time-limited pulse of the incoming photon (2t = 40fs) and 
various numbers of the layers of an absorbing barrier: AT =11 (dotted-dashed line), 
N = 41 (full line), iV = 49 (dashed line). The data of the absorbing barrier are the 
same as in Fig. |]. 
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